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Background: The SARS-CoV-2 pandemic currently prevails worldwide. To understand the immunological sig- 
nature of SARS-CoV-2 infections and aid the search and evaluation of new treatment modalities and vaccines, 
comprehensive characterization of adaptive immune responses towards SARS-CoV-2 is needed. 
Methods: We included 203 recovered SARS-CoV-2 infected patients in Denmark between April 3" and July 
9" 2020, at least 14 days after COVID-19 symptom recovery. The participants had experienced a range of dis- 
ease severities from asymptomatic to severe. We collected plasma, serum and PBMC’s for analysis of SARS- 
CoV-2 specific antibody response by Meso Scale analysis including other coronavirus strains, ACE2 competi- 
tion, IgA ELISA, pseudovirus neutralization capacity, and dextramer flow cytometry analysis of CD8* T cells. 
The immunological outcomes were compared amongst severity groups within the cohort, and 10 pre-pan- 
demic SARS-CoV-2 negative controls. 
Findings: We report broad serological profiles within the cohort, detecting antibody binding to other human 
coronaviruses. 202(>99%) participants had SARS-CoV-2 specific antibodies, with SARS-CoV-2 neutralization 
and spike-ACE2 receptor interaction blocking observed in 193(95%) individuals. A significant positive corre- 
lation (r=0.7804) between spike-ACE2 blocking antibody titers and neutralization potency was observed. 
Further, SARS-CoV-2 specific CD8* T-cell responses were clear and quantifiable in 95 of 106(90%) HLA-A2* 
individuals. 
Interpretation: The viral surface spike protein was identified as the dominant target for both neutralizing 
antibodies and CD8* T-cell responses. Overall, the majority of patients had robust adaptive immune 
responses, regardless of their disease severity. 
Funding: This study was supported by the Danish Ministry for Research and Education (grant# 0238-00001B) 
and The Danish Innovation Fund (grant# 0208-00018B) 
© 2021 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license 
(http://creativecommons.org/licenses/by-nc-nd/4.0/) 





1. Introduction 


health crisis, which is extensively affecting almost all aspects of the 
global human society. An important aspect of SARS-CoV-2 replication 


The year of 2020 has been thoroughly marked by the outbreak of 
severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) [1]. 
Reported in China December 2019, the outbreak was formally 
declared a pandemic by the WHO in March 2020 [2]. With millions of 
cases confirmed across 200 countries, the virus has claimed more 
than 2.4 million lives [3]. The SARS-CoV-2 epidemic is an ongoing 
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is binding and infection of the host cell. The viral spike protein recep- 
tor binding domain (RBD) interacts with angiotensin-converting 
enzyme 2 (ACE2), found on the cell surface, thereby mediating viral 
infection [4,5]. Coronavirus Disease 2019 (COVID-19) symptoms 
manifest primarily as a respiratory disease, with emergent complica- 
tions of several organs in cases of severe disease [6]. While efforts are 
converging globally to develop, distribute and evaluate an effective 
vaccine [7-13], our broader basic understanding of the adaptive 
immune response towards SARS-CoV-2 is still limited. 
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Research in Context 


Evidence before this study 


SARS-CoV-2 can cause severe and deadly infections. However, 
our immunological understanding of this new viral infection is 
limited, with reports focusing largely on cohorts of specific dis- 
ease severity, ages, co-morbidities etc. Recently, several vac- 
cines have been developed to help limit transmission and 
prevent the current pandemic. To understand the protective 
properties of the immune response developed during primary 
infections and to inform further vaccine development and eval- 
uation, a basic understanding of the adaptive immune response 
developed during primary SARS-CoV-2 infections is needed. 


Added value of this study 


In this context, we investigated, the adaptive immune response 
developed during SARS-CoV-2 infections in 203 recovered 
patients experiencing a full spectrum of disease severity, from 
asymptomatic infections to severe cases requiring hospitaliza- 
tion. The analysis of both binding and neutralization capacity of 
participant antibodies, alongside CD8* T-cell responses, 
towards multiple SARS-CoV-2 epitopes, provides a broad char- 
acterization of the adapted response during primary virus 
infection. We found that the vast majority of recovered individ- 
uals have clearly detectable and functional SARS-CoV-2 spike 
specific adaptive immune responses, despite diverse disease 
severities. 


Implications of all the available evidence 


The detection of both a humoral and cellular functional spike 
specific immune response in the vast majority of the individu- 
als, irrespective of asymptomatic manifestations, suggests and 
encourages further exploration of whether primary infections 
provide protection to reinfection. Furthermore, both the novel 
Meso Scale multiplex serology assay and T-cell dextramer 
staining applied here, proved powerful tools for future evalua- 
tions of vaccine efficacy. 











Several studies have described the general adaptive immune 
responses towards SARS-CoV-2, showing that SARS-CoV-2 specific B 
and T cells are generated during infections. First immunoglobulin (Ig) 
M and later IgG SARS-CoV-2 spike specific antibodies are readily 
detected in COVID-19 patients [14-17]. Evaluations by neutralization 
assays have confirmed the ability of the generated antibodies to pre- 
vent viral infections in vitro [18-20]. The limited number of confirmed 
cases suffering reinfections post recovery [21-24], and high degree of 
protective immunity against viral re-challenge shown in vivo in 
macaque challenge studies [25], suggest that the immunological 
response developed during primary infections provide at least some 
protection against reinfection. Additionally, SARS-CoV-2 specific T- 
cell activation has also been documented in a range of studies [26- 
28]. However, many studies are limited to specific disease severity 
populations, and small or none RT-PCR verified cohorts. 

Currently, in depth characterization of the adaptive immune 
response to SARS-CoV-2 in large cohorts representing the full disease 
spectrum, as well as the development of functional, and easily scalable, 
serological assays, are needed to guide and support rapid vaccine devel- 
opment and efficacy evaluation. Here, we have delineated the humoral 
and cellular immune responses in 203, RT-PCR verified, recovered 
SARS-CoV-2 patients. We evaluated the quantity and potency of anti- 
bodies in each individual towards several different coronaviruses and 
antigens, using both a SARS-CoV-2 spike pseudovirus neutralization 
assay and a novel Meso Scale Diagnostics (MSD) multiplex platform 


[29]. We further quantified the breadth and magnitude of single-epitope 
SARS-CoV-2 specific CD8* T cells, using dextramer flow cytometry. Thus, 
we report an extensive panel of adaptive immune parameters in the 
context of disease severity, to provide an outline of the general broad 
and functional SARS-CoV-2 specific adaptive immune response 
observed across the full COVID-19 disease spectrum. 


2. Methods 
2.1. Study design and sample collection 


Samples were collected from a cohort of 203 individuals who 
had recovered from COVID-19. Participants were enrolled at 
Department of Infectious Diseases at Aarhus University Hospital, 
Denmark from April 3°? to July 9" 2020. Inclusion criteria were 
as follows; 1) Age above 18 years; 2) PCR verified SARS-CoV-2 
within the preceding 12 weeks; 3) Full recovery from acute 
COVID-19 illness; 4) Able to give informed consent. Exclusion cri- 
teria were; 1) Ongoing febrile illness; 2) Immunosuppressive 
treatment and/or known immunodeficiency; 3) Pregnancy. 301 
individuals were invited to participate in the study, of which 203 
responded. All 203 responders met the study criteria and were 
included. Samples were collected at least 14 days after recovery 
and a maximum of 12 weeks after SARS-CoV-2 PCR-verified diag- 
nosis. One patient ID116 only had serum collected, and thus is 
absent from plasma neutralization and T-cell analyses. 

Individuals were allocated to three groups according to the severity 
of COVID-19 illness, based on the criteria: 1) Home/outpatient, not 
experiencing any limitations in daily activities; 2) Home/outpatient, 
certain limitations in daily activity level (fever, bedridden during ill- 
ness); 3) All hospitalized patients, regardless of need for supplemental 
oxygen treatment, or ICU admission with/without mechanical ventila- 
tion. Additional data regarding demographic and clinical characteris- 
tics of this cohort has been reported elsewhere [30]. 


2.2. Ethics 


Each participant provided informed written consent prior to any 
study activities. The study was approved by The National Health 
Ethics Committee (#1-10-72-76-20) and the Danish Data Protection 
Agency (case number not applicable) 


2.3. Serology 


IgG antibodies were measured in serum samples using the MSD 
Coronavirus Plate 1 (Cat. No. N05357A-1, Meso Scale Diagnostics, 
Rockville, Maryland), a solid phase multiplex immunoassay, with 10 
pre-coated antigen spots in a 96-well format, with an electro-chemi- 
luminescence based detection system. The SARS-CoV-2 related anti- 
gens spotted were CoV-2 Spike, CoV-2 RBD, CoV-2 NTD, and CoV-2 
nucleocapsid. The remaining spots comprised antigens from other 
respiratory pathogens: Spike protein from SARS-CoV-1, MERS coro- 
navirus, and two seasonal coronaviruses OC43, HKU1. BSA served as 
a negative control, as previously described [29]. Unspecific antibody 
binding was blocked using MSD Blocker A (Cat. No. R93AA-1). 
COVID-19 patient serum samples and control samples were diluted 
1:4630 in MSD Diluent 100 (Cat. No. R50AA-3). After sample incuba- 
tion, bound IgG was detected by incubation with MSD SULFO-TAG 
Anti-Human IgG Antibody and subsequently measured on a MESO 
QuickPlex SQ 120 Reader (Cat. No. AIOAA-0) after addition of GOLD 
Read Buffer B (Cat. No. R60AM-2). 


2.4, ACE2 competition assay 


Spike and RBD targeting antibodies with the ability to compete 
with ACE2 binding were measured using the MSD Coronavirus Plate 
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1 as above. COVID-19 blocking antibody calibrator and 1:10 diluted 
patient and control serum samples were incubated after plate block- 
ing. SULFO-Tag conjugated ACE2 was added before washing, allowing 
ACE2 to compete with antibody binding to spike and RBD antigens 
immobilized on the plate. Bound ACE2 was detected as described for 
the serology assay above, and antibody concentrations were subse- 
quently calculated using the MSD Discovery Workbench software. 


2.5. ELISA 


IgA antibodies were measured using the Anti-SARS-CoV-2 IgA 
ELISA from Euroimmun [31] (Euroimmun Medizinische Labordiag- 
nostika AG, Lubeck, Germany, Cat. No. El 2606-9601 A), according to 
manufacturer’s instructions. In brief, antibodies in serum samples 
diluted 1:200 were captured by recombinant S1 domain of SARS- 
CoV-2 spike protein immobilized in microplate wells. IgA type anti- 
bodies were detected by incubation with peroxidase labelled anti- 
human IgA followed by a chromogen solution, resulting in color 
development in positive wells. Signal was read at 450 nm with refer- 
ence measurements at 650 nm, which were used for background sig- 
nal corrections. Results were analyzed relative to the ELISA kit 
calibrator, as a ratio between sample absorbance and calibrator 
absorbance. 


2.6. Cells and plasmids 


All cell lines were incubated at 37°C and 5% CO, in a humidified 
atmosphere. BHK-G43 cells, previously described and kindly pro- 
vided by M. Hoffmann and S. PohImann [32,33], were cultured in Dul- 
becco’s modified eagle’s medium (DMEM) (Biowest VWR Cat no. 
L0101-500), containing 5 % Fetal Bovine Serum (FBS) (Biowest prod- 
uct code: $181H-500) and 50 U/mL Penicillin G/Streptomycin (P/S) 
(Biowest, product code: LO022-100) where Zeocin (100jzg/ml)(Inviv- 
ogen cat.code: ant-zn-1) and Hygromycin (50jg/ml)(Corning, Fisher 
Scientific product code: 15313681) were added at every fourth pas- 
sage. Induction of VSV-G glycoprotein was performed with 1078M 
mifepristone (Merck cat.no M80-46). HEK293T cells (RPID: 
CVCL_0063) were cultured in DMEM, containing 10% FBS and 
50 U/mL P/S. Vero76 cmyc hTMPRSS2 cells kindly provided by M. 
Hoffmann and S. Pohlmann [4] were cultured in DMEM supple- 
mented with 10% FBS, 50 U/mL P/S, and 10 ug/mL Blasticidin(Gibco, 
REF: A1113903). 

The construction of pCG1-SARS-2-Spike has been previously 
described and was kindly provided by M. Hoffmann and S. Pohlmann 
[4,34]. Briefly, SARS-2-S from the Wuhan-Hu-1 isolate (NCBI Ref.Seq: 
YP_009724390.1) coding sequence was PCR-amplified and cloned 
into the pCG1 expression vector via BamHI (ThermoFisher Cat.no: 
ERO051) and Xbal (ThermoFisher Cat.no: ERO685) restriction sites. 


2.7. Virus production 


For generation of VSV*AG(luc)-G particles BHK-G43 cells were 
seeded day 1 to reach a confluence of 70-80% at day 2, where Mifep- 
ristone (1078 M) was added to induce transcription of glycoprotein G. 
After 6 hours the medium was replaced with fresh DMEM containing 
5% FBS, 50 U/mL P/S, and VSV*AG(luc) (kindly provided by M. Hoff- 
mann and S. PohIman [32]) at MOI = 0.3. After 1 hour of incubation at 
37°C BHK-G43 cells were washed three times in PBS and fresh media 
was added. Cells were incubated for 24 hours, after which the super- 
natant was centrifuged at 2000 xg for 10 min at room temperature to 
pellet cellular debris, and stored at -80°C. 

VSV*AG(luc)-SARS-2-S pseudovirus was produced by transfection 
with pCG1-SARS-2-S followed by transduction with VSV*AG(luc)-G. 
HEK293T cells were seeded in DMEM containing 10% FBS and 50 U/mL 
P/S to reach 70-80% confluence the next day. 2ug plasmid was used 
per 1 x 10° cells and incubated with PEI (3:1) for 30 min at room 


temperature. The transfection mixture was added to the cells, and 
incubated for 18 hours at 37°C. Cells were washed twice with PBS, 
transduced with VSV*(luc)+G at MOI = 2, and incubated for 2 hours. 
The virus was removed by gently washing with PBS twice, and fresh 
DMEM containing 10% FBS and 50 U/mL P/S was added. Cell superna- 
tant was harvested after 24 hours, centrifuged at 2000 xg for 10 min to 
eliminate cellular debris, and stored at -80°C immediately. A VSV*AG 
(luc)-mock was generated simultaneously to allow subtraction of any 
remaining background from VSV*AG(luc)-G signals. 


2.8. Neutralization assay 


The SARS-CoV-2 neutralization capacity of plasma was assessed 
through infection of Vero76 cmyc hTMPRSS2 cells, with VSV*AG(luc)- 
SARS-2-S pseudovirus particles. Neutralization was conducted as fol- 
lows: Plasma samples were thawed and heat-inactivated at 56°C for 
45 min. Subsequently, five-fold serial dilutions of each plasma sample 
were made in DMEM containing 10% FBS and 50 U/mL P/S. 25 uL of 
each dilution for all plasma samples was incubated with 50 uL 
VSV*AG(luc)-SARS-2-S at MOI = 0.01 for 1 hour at 37°C, in a flat bot- 
tomed 96-well plate. Successively, 20,000 Vero76 cmyc hTMPRSS2 
cells, in 50 uL DMEM containing 10% FBS and 50 U/mL P/S were 
added to each well, and incubated at 37°C for 20 hours. The final total 
plasma dilutions tested ranged from 1:25 — 1: 1953125, and were all 
analyzed in duplicates. Cells were prepared for flow cytometry by 
first gently removing the culture media, and washing once with PBS. 
Second, cell suspensions were made by incubating each well with 
75uL Trypsin with EDTA(Biowest product code: X0930-100) for 
15 min at 37°C, followed by centrifugation at 500 g for 5 min at room 
temperature, and re-suspension in DMEM containing 10% FBS and 
50 U/mL P/S. Cells were fixed in 1% PFA(VWR Cat. No: AA43368-9M) 
for at least 15 min at 4°C, before eGFP expression was analyzed using 
a Miltenyi Biotec MACSquant16 flow cytometer. The VSV*AG(luc)- 
mock eGFP background signal was subtracted from all samples. 


2.9. HLA-A2 typing and dextramer staining by flow cytometry 


For HLA-A2 typing 2 million cryopreserved PBMCs were thawed, 
stained at room temperature for 20 min with HLA-A2 (clone BB7.2, 
Biolegend Cat. No. 343328) or matching isotype control (Biolegend 
Cat. No. 400356) and acquired on a five-laser Fortessa flow cytome- 
ter. The dextramer stains were then performed on the HLA-A2 posi- 
tive samples as follows. PBMCs were incubated at room temperature 
for 30 min with the following SARS-CoV-2 dextramers (all from 
Immundex): A*0201/TLACFVLAAV-PE (Cat. No. WB3848-PE), A*0201/ 
GMSRIGMEV-FITC (Cat. No. WB5751-FITC), A*0201/LLLDRLNQL-APC 
(Cat. No. WB5762-APC), A*0201/ILLNKHIDA-PE (Cat. No. WB5848- 
PE), A*0201/RLNEVAKNL-FITC (Cat. No. WB5750-FITC), A*0201/ 
YLQPRTFLL-APC (Cat. No. WB5824-APC), A*0201/VLNDILSRL-PE (Cat. 
No. WB5823-PE), A*0201/NLNESLIDL-FITC (Cat. No. WB5850-FITC), 
A*0201/FIAGLIAIV-APC (Cat. No. WB5825-APC), A*0201/LLLNCLWSV- 
PE (Cat. No. WB3513-PE), or positive/negative control dextramers: 
A*0201/NLVPMVATV-PE (Cat. No. WB2132-PE, Pos. Control, CMV), 
A*0201/NLVPMVATV-FITC (Cat. No. WB2132-FITC, Pos. Control, 
CMV), A*0201/NLVPMVATV-APC (Cat. No. WB2132-APC, Pos. Control, 
CMV), A*0201/Neg. Control-PE (Cat. No. WB2666-PE), A*0201/Neg. 
Control-FITC (Cat. No. WB2666-FITC), A*0201/Neg. Control-APC (Cat. 
No. WB2666-APC). Cells were washed and stained with viability dye 
(Zombie Violet, Biolegend, Cat. No. 423114) and CD8 (Clone RPA-T8, 
BD, Cat. No. 563795) and acquired on a five-laser Fortessa flow 
cytometer. 


2.10. Data and statistical analyses 


Flow cytometry data was analysed using FlowJo (Version 
10.7.1). All data was processed and graphed in GraphPad Prism 
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version 8.4.3. For the Meso Scale serology and ACE2 competition 
assays, and ELISA IgA assay, the cohorts chemiluminescence sig- 
nal or antibody titer was compared to the same variable for 10 
negative controls, within each antigen target separately, using a 
Mann-Whitney test. Mann-Whitney tests were also used to com- 
pare variables between the cohorts severity groups, comparing 
each severity group with each other individually, in the cases of: 
chemiluminescent signal towards SARS-CoV-2 spike and nucleo- 
capsid, SARS-CoV-2 spike pseudovirus neutralization IC50 values, 
ACE2 blocking SARS-CoV-2 antibodies in ng/ml, and the cumula- 
tive CD8 T-cell response %. Spearman’s rank correlation analysis 
was used to access the correlation between SARS-CoV-2 spike 
pseudovirus neutralization IC50 values and ACE2 blocking SARS- 
CoV-2 antibodies in ng/ml, as well as correlations between the 
cumulative CD8 T-cell response % towards SARS-CoV-2 A2-Pepti- 
des, and IC50 values as well as Meso Scale serology chemilumi- 
nescence signals and ACE2 blocking antibody titers towards the 
epitopes listed in table 2. Neutralization curves were plotted with 
three parameter non-linear fits, from which IC50 values were cal- 
culated. All values in the heatmap was normalized to values 
between 1 and 0, based on the highest and lowest value for each 
individual variable displayed. No randomization or blinding was 
carried out, as patient enrolment was based on a history of 
COVID-19 diagnosis, with disease severity being unknown until 
enrolment. Due to the number of comparisons made, multiple 
testing was not accounted for. No statistical corrections were 
made to account for potential confounding. p < 0.05 was inter- 
preted as statistically significant. P-values are indicated as fol- 
lows: n.s. = not significant,* = p < 0.05, ** = p < 0.01, *** = p < 
0.001, **** = p < 0.0001. 


2.11. Role of the funding source 


This study was supported by a grant from the Danish Ministry for 
Research and Education (grant# 0238-00001B) and The Danish Inno- 
vation Fund (grant# 0208-00018B). None of the funding sources had 
any role in any aspect pertinent to the study. The corresponding 
author had access to all data in the study and held the final responsi- 
bility for the decision to submit for publication. 


3. Results 
3.1. Patient enrollment 


We studied the adaptive immune response towards SARS-CoV-2 
among 203 Danish patients who had recovered from COVID-19. We 
have recently described the cohorts clinical characteristics thor- 
oughly [30] a basic overview of which is shown in Table 1. The 
median age of individuals was 47 years (range: 21 — 79), and 45% 


Table 1 
Cohort characteristics. 





Demographics and Clinical Characteristics at Baseline 








Characteristics n=203 

Age, years, median (range) 47 21-79) 

Female sex, no (%) 92 45) 

HLA-A2*, no (%) 113 56) 

COVID-19 disease severity, no (%) 

Home/outpatient, no limitation of daily activities 17 8) 
(asymptomatic/mild) 

Home/outpatient, limitation of daily activities 152 75) 
(moderate) 

Hospitalized (severe) 34 17) 

Duration of COVID-19 symptoms, days, median 13 0-68) 
(range) 

Time from recovery to inclusion, day, median (range) 31 14-61) 





were female. The cohort was divided into three COVID-19 disease 
severity groups. 1: Home/outpatients with no limitation of daily 
activities (8%), 2: Home/outpatients with a limitation of daily activi- 
ties (75%), and 3: Hospitalized patients (17%). The median duration of 
COVID-19 symptoms was 13 days (range: 0 — 68). Enrolment 
occurred at least 14 days after the end of COVID-19 related symp- 
toms, with a median of 31 (range: 14 — 61) days from time of recov- 
ery to study enrolment. To allow comparison of immunological 
outcomes from SARS-CoV-2 infection recovered patients, samples 
from 10 healthy Danish individuals enrolled in a study conducted 
prior to the current COVID-19 pandemic were included as controls 
[35]. 

All individuals were assigned a COVID-19 severity group depend- 
ing on their course of disease. Group 1 consisted of asymptomatic 
individuals with no limitations in their daily activities. Group 2 of 
moderately sick, able to recover at home. Finally, group 3 comprises 
all hospitalized individuals, including those with/without oxygen 
requirement and/or ICU admission. 


3.2. Human coronavirus serology 


First, we analysed the presence of IgG antibodies towards multiple 
human coronaviruses in serum, using the multiplex MSD platform. 
Compared to controls, we found significantly elevated levels of IgG 
antibodies in spike RBD, spike N-terminal domain (NTD), and the 
nucleocapsid (p<0.0001, Fig. 1a). Furthermore, IgG antibodies from 
SARS-CoV-2 infected individuals exhibited strongly increased reac- 
tivity towards spike protein from other human beta coronaviruses: 
SARS-CoV-1 and Middle East respiratory syndrome (MERS), as com- 
pared to the negative controls. Further, increased IgG levels towards 
the seasonal beta coronavirus strains: HKU1 and OC43, compared to 
IgG from the control group were also observed (p<0.0001, Fig. 1b). 
No difference was detected in IgG levels to the negative bovine serum 
albumin (BSA) control between SARS-CoV-2 patients and controls. 
Importantly, 202 out of the 203 individuals analysed here, developed 
detectable antibodies, otherwise absent in the historical controls, 
against both full-length SARS-CoV-2 Spike and RBD antigens, during 
SARS-CoV-2 infections. Likewise, robust production of IgA antibodies 
was also observed for nearly all infected individuals, with SARS-CoV- 
2 spike specific IgA levels being significantly elevated compared to 
controls in 201 of the 203 individuals (Fig. 1c). Additionally, SARS- 
CoV-2 IgG levels towards both spike and nucleocapsid antigens, cor- 
related positively with the disease severity. (Fig. 1d+e). Overall, we 
conclude that more than 202(99%) of the SARS-CoV-2 infected indi- 
viduals in this cohort had readily detectable antibodies to SARS-CoV- 
2 spike antigen, and that broad IgG immunological recognition of 
SARS-CoV-2 as well as several different coronavirus develops during 
COVID-19. Additionally, the magnitude of spike-targeting antibodies 
increases with disease severity. 


3.3. SARS-CoV-2 pseudovirus neutralization 


Next, we investigated the functional neutralization capacity of total 
plasma antibodies in vitro, using VSV pseudotyped virus expressing 
SARS-CoV-2 spike protein. Antibody neutralizing potency was evalu- 
ated by applying serial dilutions of plasma, yielding infectivity titration 
curves for each of the SARS-CoV-2 infected individuals and the controls 
(Fig. 2a). We found that 95.5% of the individuals (193 of 202) were able 
to neutralize SARS-CoV-2 spike pseudoviruses with 100% inhibition at 
the lowest (1:25) plasma dilution. IC50 values were calculated from the 
neutralization curves, and assigned to each individual as a measure of 
antibody neutralization potency. Serum from the remaining nine indi- 
viduals (4.5%) was unable to fully neutralize viral infection, producing 
neutralization curves comparable to that of the uninfected controls. No 
legitimate IC50 value could be calculated for these individuals, and con- 
sequently they were excluded from further analyses using this 


S.S. Nielsen et al. / EBioMedicine 68 (2021) 103410 5 


SARS-CoV-2 IgG 


hk kkk paee pose 





10° Fa 2 


i 
104 





z 
+Ë, s 


F 


Electro Chemiluminescence 
è 
i 





Patients Controls Patients Controls Patients Controls Patients Controls 





Full Length Spike Spike - RBD Spike -NTD Nucleocapsid 


Human Coronavirus Spike IgG 


kikk tik tek 


kee 


SARS-CoV-2 Spike IgA 
102 


kk 


kkk 


=à 
o 





IgA Ratio 
S 


è 
SARS-CoV-2 Spike specific 
Electro Chemiluminescence 





Patients Controls 1 2 3 


Severity Group 


kkk 


kkrk 





Electro Chemiluminescence 





10° 


105 





104 


SARS-CoV-2 Nucleocapsid specific 
Electro Chemiluminescence 





Patients Controls Patients Controls Patients Controls Patients Controls Patients Controls 1 2 3 


Ll 
SARS-CoV-1 





MERS 


Severity Group 


Fig. 1. Extensive IgG and IgA presence with multiple SARS-CoV-2 antigens. a+b) Serum IgG levels for all individuals (n=203) and 10 pre-pandemic healthy controls. IgG was detected 
against SARS-CoV-2 Spike, RBD (receptor binding domain), NTD (N-terminal domain), nucleocapsid and non-SARS-CoV-2 spike proteins of other corona virus. Data are blank-cor- 
rected electro chemiluminescent signal measured by MSD multiplex serology assays. c) Serum IgA levels for all individuals (n=203) and eight pre-pandemic healthy controls, mea- 
sured by ELISA. IgA is shown as a ratio against a standard calibrator. d+e) Distribution of IgG volumes between each disease severity group, for both SARS-CoV-2 spike (d) 
* = p=0.0416 and nucleocapsid (e) * = p=0.0271. Data are blank-corrected electro chemiluminescent signal measured by MSD multiplex serology assays. Scatter plots with individual 
data points are shown with median (wide line) and interquartile range (narrow lines). Statistical comparison between groups were done by Mann-Whitney test. n.s = not significant, 


=" = p<0.0001. 


parameter. Collectively, the IC50 values of all 193 neutralizing individu- 
als span evenly across four orders of magnitude (Fig. 2b). In concur- 
rence with the analysis in Fig. 1d+e, we observed lower IC50 values 
among individuals experiencing mild symptoms compared to those 
with moderate (p=0.0008) or severe COVID-19 (p<0.0001) (Fig. 2c). 
We conclude that in this large cohort, with considerable diversity in 
disease severity, the vast majority (>95 %) of SARS-CoV-2 infections 
lead to the production of effective neutralizing antibodies, and that 
neutralization potency increases with disease severity. 


3.4. Antibodies efficiently block ACE2 receptor binding 


We continued the characterization of SARS-CoV-2 antibody func- 
tionality, using an MSD SARS-CoV Spike — ACE2 competition assay 
(Fig. 3a). This allowed us to measure the quantity of antibodies able 
to block the interaction between the ACE2 receptor and SARS-CoV-2 
full-length spike protein, SARS-CoV-2 RBD, and SARS-CoV-1 spike 
protein. Many of the recovered individuals reached the assay’s upper 
limit of quantification, and a clear increase in the quantities of serum 
ACE2 blocking antibodies was observed for all three antigens com- 
pared to historic controls (p<0.0001) (Fig. 3b). The levels of antibod- 
ies blocking SARS-CoV-2 Spike — ACE2 receptor interaction was 
increased in >99% of the individuals (202 of 203) compared to unin- 
fected controls. The individual antibody concentrations also 


correlated to the time from disease recovery to inclusion (S1 Fig. 1b). 
Nevertheless, we found that those experiencing severe COVID-19 
had significantly greater levels of SARS-CoV-2 spike specific ACE2 
blocking antibodies, compared to individuals with mild to moderate 
disease (p<0.0001, Fig. 3c). Both the pseudovirus cell-based neutrali- 
zation assay and the SARS-CoV Spike — ACE2 competition assay 
investigate the presence of functional antibodies towards SARS-CoV- 
2. We identified a highly significant correlation between the IC50 val- 
ues from the pseudovirus neutralization assay and the concentration 
of SARS-CoV-2 spike specific antibodies capable of blocking ACE2 
receptor interaction (p>0.0001 Fig. 3d). In conclusion, we observed 
that nearly all individuals produce antibodies that target the spike 
protein-ACE2-receptor interaction and that the level of these anti- 
bodies was increased with severe disease. Further, the virus neutrali- 
zation capacity increased in conjunction with the amount of 
functional ACE2 blocking antibody present in serum. 


3.5. Collected serological analysis 


Next, we constructed a heatmap compiling all humoral immuno- 
logical data, to gain a cohort wide perspective of the overall antibody 
response developed during SARS-CoV-2 infection. We ranked indi- 
viduals according to their antibody response potency from the pseu- 
dovirus neutralization assay (IC50 value), displaying their respective 
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Fig. 2. SARS-CoV-2 neutralization capacity correlates with disease severity. a) Representative neutralization curves for control ID308, and individuals ID54, ID194, and ID203, quan- 
tified as eGFP* cells by flow cytometry. Control plasma was unable to neutralize below a 50% infection rate, where SARS-CoV-2 recovered patients accomplish 100% neutralization 
at the lowest plasma dilution. X-axis shows the log10 transformed patient plasma dilution, from 1:25 — 1:1,953,125. Error bars represent mean and s.e.m. of duplicate determina- 


tions. Three-parameter non-linear fit is plotted. b) IC50 values calculated from neutralization curves, graphed from lowest (left) 


— highest (right) within the cohort. Error bars show 


95% confidence interval. Nine individuals unable to neutralize 100% are represented with the value zero on the y-axis far left, n = 202. c) Distribution of IC50 values between disease 
severities. *** = p=0.0008. Scatter plot with individual data points shown with median (wide line) and interquartile range (narrow lines). Statistical comparison were by Mann-Whit- 


ney test. **** = p < 0.0001, n = 193. 


immunological variables underneath (Fig. 4). We observed, that the 
neutralization capacity was clearly linked to the overall antibody lev- 
els present in the patients. Interestingly, it was further evident, that 
the best (top 10%) neutralizers of the cohort displayed a correspond- 
ing increase in the overall breadth of their antibody response, 
towards all the investigated coronavirus antigens. Importantly, 
strong pseudovirus neutralization profiles were almost exclusively 
seen in individuals with antibodies that potently block spike-ACE2 
receptor interaction. We therefore conclude that the best neutralizers 
exhibit a broader variety of antibodies and have greater levels of 
spike binding and receptor-blocking antibodies. 


3.6. Epitope specific CD8+ T cell-responses 


We then went on to explore the epitope specific T-cell responses 
in SARS-CoV-2 recovered individuals. We analysed the reactivity of 
CD8* T cells from 106 HLA-A2* individuals in the cohort for their 
specificity to nine different SARS-CoV-2 epitopes using dextramer 
staining flow cytometry (Fig. 5a). Overall, Membraneg,-7o (M) (epi- 
tope 1), Nucleocapsidz22-230 (N) (epitope 3), and Spike269-277 (S) (epi- 
tope 6) were the most commonly recognized epitopes with positive 
responses detected in 19(17%), 27(25%) and 86(81%) of individuals, 
respectively (Fig. 5b). Interestingly, these three epitopes originate 
from three separate SARS-CoV-2 proteins (Fig. 5a). The frequency of 
SARS-CoV-2 specific CD8* T-cells was similar across all nine HLA-A2* 
epitopes tested, with the highest individual responses observed for 
N222-230 and S269-277 (epitopes 3 and 6) (Fig. 5c). Only 10% of the HLA- 
A2* individuals (11 of 106) had no detectable response to any of the 
epitopes tested, while the remaining 90% responded to at least one, 
and up to seven, of the analysed epitopes. (Fig. 5d). We compared the 
cumulative frequency of SARS-CoV-2 specific CD8* T cells across the 
disease severity groups and observed no significant difference 
(Fig. 5e). However, we did observe significant albeit weak correla- 
tions between the cumulative frequency of SARS-CoV-2 specific CD8* 


T-cells and the majority of the serological immunological parameters 
analysed, including pseudovirus neutralization IC50 values as well as 
SARS-CoV-2 specific antibody production and ACE2 blocking ability, 
as outlined with correlation coefficients in table 2. Additionally, the 
individual cumulative frequency of SARS-CoV-2 specific CD8* T-cells 
also correlated to the time from disease recovery to inclusion. (S1 
Fig. 1a) The 11 individuals with no detectable CD8* T-cell responses 
were evenly distributed among the disease severity groups and dis- 
played varying antibody neutralization capacity (S1 Fig. 3). Based on 
this we were only able to identify two individuals with both no 
detectable neutralizing antibodies and no detectable CD8* T-cell 
responses. Thus, we conclude that 90% of SARS-CoV-2 infected indi- 
viduals mount a detectable CD8* T cell response, towards the nine 
epitopes tested, irrespectively of disease severity. We further con- 
clude that the broadest targeted epitope in this cohort is located in 
the spike protein. Lastly, there is an overall weak but statistically sig- 
nificant correlation of antibody responses and CD8* T-cell responses. 

Spearman’s rank coefficient correlations displaying the relation- 
ship between the overall magnitude of CD8* T-cell responses toward 
SARS-CoV-2 epitopes in %, and neutralization IC50 values, IgG chemi- 
luminescent signal and ACE2 blocking antibody titers in ng/ml, for all 
SARS-CoV-2 antigens investigated. 


4. Discussion 


We aimed to characterize the cellular and humoral adaptive 
immune response in a large cohort of RT-PCR verified SARS-CoV-2 
recovered patients, spanning a full spectrum of COVID-19 severity. 
Overall, our results show that the majority of patients developed a 
robust and broad both humoral and cellular immune response to 
SARS-CoV-2. However, our data may also help explain that some rare 
individuals have no detectable immunological memory to SARS-CoV- 
2, and may therefore be at risk of re-infection as it has been reported 
in a few case reports [21-24]. 
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Cumulative CD8* T-cell responses in correlation to serology. 





Correlations to cumulative epitope specific CD8+ T-cell responses 





Immunological parameter r-value p-value 
1C50 values 0.2542 0.0107 
SARS-CoV-2 Spike ACE2 Blocking antibodies ng/mL 0.2906 0.0147 
SARS-CoV-2 RBD ACE2 Blocking antibodies ng/mL 0.3057 0.0101 
SARS-CoV-2 Spike IgG 0.2659 0.0261 
SARS-CoV-2 RBD IgG 0.2704 0.0236 
SARS-CoV-2 N-Terminal Domain IgG 0.2918 0.0143 
SARS-CoV-2 Nucleocapsid IgG 0.2102 0.0807 


confirmed the functionally neutralizing and ACE2 blocking capabili- 
ties of the SARS-CoV-2 spike and RBD specific antibodies. Noticeably, 
this infers the development of a robust humoral immune response 
within the vast majority of the COVID-19 recovered population. Fur- 
thermore, nearly all individuals also have SARS-CoV-2 specific IgA 
responses, clearly indicating functional rigorous class switching and 
antibody maturation. This presence of IgA is crucial for the immuno- 
logical protection at mucosal barriers, and hence protection against 
future SARS-CoV-2 exposures. 

All serological and functional data collected show that both anti- 
body levels and neutralization potency correlate significantly with 
the disease severity. This indicates that severe disease manifestation 
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is not caused by a lack of adaptive immunity, which is in line with 
previous reports [37,38]. Hence, we suggest that the prolonged dis- 
ease course, and consequent larger exposure to virus experienced in 
hospitalized patients, may provide a timeframe in which enhanced 
antibody affinity maturation takes place, compared to shorter course 
mild infections. Of note the duration of illness, also impact the time 
from infection to sampling in this study as participants had to be 
completely free of symptoms prior to hospital visit. This time from 
infection detection to sample collection is a clear limitation of this 
study, because disease duration is heterogeneous within the cohort. 
In the setting of the first-wave of the pandemic, we could not allow 
individuals with e.g. ongoing cough to enter the hospital facilities. As 
COVID-19 has a very broad disease spectrum duration, some partici- 
pants had to wait longer for their visit until the lingering acute symp- 
toms had resolved completely. Alternatively, we would not have 
been able to report on a recovered cohort 

Studies are conflicted on the degree to which cross-reactive 
immunity between different coronavirus develop during SARS-CoV-2 
infections [18,37,39-41]. The considerable diversity of antigen recog- 
nition independent of COVID-19 severity shown here, demonstrates 
that at least some immunological cross-recognition of several differ- 
ent coronavirus is developed during SARS-CoV-2 infections. This is in 
line with data on cross-reactivity in CD4* T-cell epitopes between 
seasonal coronaviruses and SARS-CoV-2 [42]. The cross-reactivity 
observed between SARS-CoV-2, SARS-CoV-1 and MERS, may be due 
to conserved epitopes between these viruses, as prior infections with 
SARS-CoV-1 or MERS within our Denmark based cohort are highly 
unlikely. Such potential cross-reactivity could arise through either 
newly generated SARS-CoV-2 specific antibodies reacting with con- 
served epitopes, or by reactivation of memory cells originally gener- 
ated against seasonal coronaviruses. The increased antibody 
recognition between the cohort and the controls for the seasonal 
coronavirus strains HKU1 and OC43, while statistically significant, 
was much lower than that observed for the remaining coronavirus 
epitopes. This suggests that the SARS-CoV-2 negative controls may 
have prior exposure to these strains, increasing the overall back- 
ground for these epitopes significantly. Importantly, the multiplex 
serological analyses we performed do not provide insight into the 
SARS-CoV-2 antibody response on a monoclonal antibody level. Here, 
further studies are needed to determine the possible protective and 
cross-reactive properties of single-antibody specificities. 

We functionally verified the antibody responses in all individuals, 
using two separate assays. The cell-based neutralization assays are at 
present the standard method for determining SARS-CoV-2 neutraliz- 
ing antibody potency. We additionally used the MSD novel coronavi- 
rus multiplex assay, recently reported by Johnson et al [29] to 
determine the ACE2 blocking capability of individual serum antibod- 
ies. The significant correlation between the two assay readouts iden- 
tifies the plate format ACE2 competition assay as a powerful, high- 
throughput, screening tool, with applications in both SARS-CoV-2 
therapeutic neutralizing antibody development, and assessments of 
functional protective antibody induction post vaccination. An 
immense global effort is currently undertaken to develop, distribute 
and evaluate effective vaccines against SARS-CoV-2, the majority of 
which are focused on inducing spike or RBD antigen specific immu- 
nity [43]. Here we demonstrate that SARS-CoV-2 spike specific, ACE2 
blocking antibodies are found in the majority of infected individuals. 
Their extensive induction, even in short-term, asymptomatic infec- 
tions, align with current vaccines designs inducing protective immu- 
nity based on spike antigens [8-13,44,45]. Nevertheless, the 
protective effect of antibodies elicited during natural infections, 
remains to be determined. 

We further report, with single-epitope resolution, a SARS-CoV-2 
specific CD8* T-cell response in 90% of the HLA-A2* individuals ana- 
lysed. This corresponds well with other studies reporting CD8* T cell 
activation in 70%—100% of recovered patients using full protein 


overlapping peptide stimulation [26,46]. The location of the top three 
immunogenic epitopes within separate proteins in the viral prote- 
ome additionally reinforces our conclusion that a broad immune 
response is generated towards SARS-CoV-2 in the general infected 
population. T-cell immunity to SARS-CoV-1 is known to persist for up 
to six years [47,48], underlining the importance of developing protec- 
tive cell based immunity to SARS-CoV-2 if long term viral protection 
is to be obtained. As an important point, the most broadly recognized 
CD8* T-cell epitope (S2¢9-277) within our cohort (responses detected 
in 81% of HLA-A2* individuals) is located in the spike antigen. Thus, 
such epitope specificity can clearly be used to assess CD8* T-cell 
immunity in the evaluation of the vaccination efforts currently 
underway. 

Surprisingly, we found that the cumulative CD8* T-cell 
response, across all epitopes, did not vary by disease severity in 
contrast to what we, and others [49], observed with antibody lev- 
els. While the limited coverage of epitopes investigated here may 
influence this observation, recent evidence suggests that persis- 
tent viral replication in otherwise recovered patients may be 
linked to CD8* T-cell response magnitude [30]. Despite the differ- 
ent observations with regard to immune responses and disease 
severity, we found overall significant relationships between 
humoral and T-cell based immunity, but all of modest strength. A 
possible explanation could be the synchronized waning of the 
magnitude of response for both immune parameters during the 
time from recovery to study enrolment. 

Of note, the use of dextramer staining is limited by inclusion of 
selected epitopes and HLA type only, and conclusions are conse- 
quently limited to the relative low epitope coverage. However, the 
advantages of the dextramer technology are superior sensitivity and 
a high degree of specificity. In the light of the relative low proteome 
coverage, the fact that only 10% of the investigated individuals did 
not have a detectable CD8* T-cell response clearly indicate a strong 
cytotoxic T-cell component in the immune response towards SARS- 
CoV-2. Furthermore, as our observations of breadth and magnitude 
in relation to the distribution of distinct SARS-CoV-2 antigens are 
similar to others [46,49] we conclude that the panel of dextramers 
applied here provide a new and sensitive representation of the gen- 
eral CD8* T-cell response to SARS-CoV-2 that will be an important 
tool in assessing long-term immunity following primary infection or 
vaccination. 

In conclusion, we observed that disease severity is closely related 
to the potency and breadth of the antibody response towards SARS- 
CoV-2. Furthermore, we identified the SARS-CoV-2 spike protein as a 
target of adaptive immunity in >99% of the cohort, irrespective of 
COVID-19 symptom manifestation. Only two individuals (<2%) had 
neither antibodies with virus neutralization capacity, nor detectable 
CD8* T-cell responses. Hence, we conclude that regardless of COVID- 
19 severity, a robust adaptive immune response towards SARS-CoV-2 
is elicited during primary infections. 
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